Memory Imperfections in Atomic Ensemble-based Quantum Repeaters 
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Quantum repeaters promise to deliver long-distance entanglement overcoming noise and loss in 
realistic quantum channels. A promising class of repeaters, based on atomic ensemble quantum 
memories and linear optics, follow the proposal by Duan et al [Nature 414, 413, 2001]. Here we 
analyse this protocol in terms of a very general model for the quantum memories employed. We 
derive analytical expressions for scaling of entanglement with memory imperfections, dark counts, 
loss and distance, and apply our results to two specific quantum memory protocols. Our methods 
apply to any quantum memory with an interaction Hamiltonian at most quadratic in the mode 
operators and are in principle extendible to more recent modifications of the original DLCZ proposal. 

PACS numbers: 03.67.Ac, 03.67.Hk, 42.50.Gy 



I. INTRODUCTION 

If one attempts to transfer quantum information by di- 
rect transmission, the communication rate decreases ex- 
ponentially with distance, due to decoherence and loss. 
Quantum repeaters achieve subexponential scaling by 
generating entanglement locally in parallel for many seg- 
ments of a short length Lq and subsequently extend- 
ing the distance by entanglement swapping until the full 
channel length L = 2 n Lo is reached d- Quantum memo- 
ries play a crucial role in quantum repeaters because op- 
erations take place on many segments in parallel and the 
operations in each segment may fail with a large probabil- 
ity. It is therefore essential to have a quantum memory 
where entanglement successfully generated in one seg- 
ment may be stored, while entanglement generation and 
connection is being attempted in other segments. One 
promising repeater protocol based on storage of light in 
atomic ensembles was proposed by Duan, Lukin, Cirac 
and Zoller (DLCZ1 [2]| and later improved upon in a 
number of papers 0,13, H, S, 0|- Extensive experimen- 
tal progress have been made toward the realisation of 
this protocol @, ©, E3, El ■ 

In this paper, we consider the effect of memory im- 
perfections in a repeater architecture closely resembling 
the original DLCZ proposal. Using a general model for 
memories, we investigate how the repeater performance 
depends on the memory properties in the presence of re- 
alistic errors, i.e. lossy fibres, detection inefficiency and 
dark counts. We then apply the results to specific en- 
semble based memories 0, E2, E3 and evaluate the per- 
formance of quantum repeaters based on these memo- 
ries. The methods we develop, though applied here to 
the DLCZ architecture, could be extended also to the 
more recent protocols (& 0, IE IE 0, El • 



II. MODELLING THE REPEATER 

The repeater that we shall consider is defined by the 
setups for entanglement generation and connection, illus- 
trated in Fig. [TJ Although the use of repeaters is moti- 



vated by the presence of errors in transmission, it is easier 
to understand the basics of the protocol in the absence 
of errors. Hence we consider this ideal case first. 

To generate entanglement in one segment of the re- 
peater, two non-degenerate parametric down converters 
(PDC's) are used. The two-mode squeezed state gener- 
ated by a single PDC with small squeezing parameter r 
is 
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fc=0 

where jfcijfe) denotes a Fock state with photons num- 
bers fci and &2 in modes 1 and 2. A single excitation is 
subtracted non-locally from the two squeezed pairs, by 
mixing one mode from each pair on a balanced beam 
splitter and conditioning on a single click (Fig. QJ,). In 
the ideal case of noiseless operations and photon number 
resolving detectors, this leaves a Bell state to be stored 
in the atomic memories 
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To extend the entanglement distance, two neighbouring 
segments are connected by mixing one mode from each 
on a balanced beam splitter and conditioning on a sin- 
gle click (Fig. [TL). After the connection, the remaining 
modes are again in a Bell state, and thus under ideal 
conditions the protocol generates a maximally entangled 
pair over the distance L. 

In Ref . [2] it was shown that the type of protocol above 
works also in the presence of photon loss. By condition- 
ing on a click one essentially purifies the state from losses. 
Other errors, however, are not purified in this way. Our 
goal here is to make a detailed investigation of the effect 
of various types of noise on the efficiency of the protocol. 
The noise sources that we will consider are: (i) transmis- 
sion losses, (ii) detector dark counts and (iii) memory 



To model the transmission losses (i), virtual beam 
splitters are inserted into the setup as illustrated in 
Fig. [TJ For simplicity we assume that the memory is 
located close to the PDC's, such that the loss in the 




FIG. 1: (a) Entanglement generation. Two PDC's (large 
boxes) emit two-mode squeezed states. One mode from each 
pair is stored in an atomic memory (circles), while the re- 
maining modes are mixed on a 50/50 beam splitter and mea- 
sured. Observation of a single click leads to entanglement 
between the atomic modes. The dashed elements are virtual 
beam splitters and PDC's modelling losses and detector dark 
counts, (b) Entanglement connection is achieved by reading 
out the connecting ends of two entangled pairs and mixing 
on a 50/50 beam splitter. A single click heralds successful 
connection. For perfect memories, this setup is equivalent to 
the DLCZ setup 0. 



memory arm is small and can be treated perturbatively 
as a memory imperfection. In the detector arms of the 
entanglement generation and connection setups, the loss 
probabilities are p gen and p con respectively. Photon loss 
due to detector inefficiency can be included in the trans- 
mission loss. 

To model the dark counts (ii) , we assume that the sig- 
nal to be measured is mixed with a thermal state, which 
we generate by inserting a virtual PDC (see Fig. [IJ . We 
choose a thermal distribution for the dark-counts, be- 
cause this is easily treated as a Bogoliubov transforma- 
tion (see below). The repeater setup presented here is 
only feasible for fide "C 1, where hd c is the average num- 
ber of dark counts per detector in one measurement cy- 
cle. Therefore only the first-order contribution from dark 
counts is considered, and the actual distribution is not 
important. The reflection coefficients of the virtual beam 
splitters are given by the photon loss probabilities, and 
the detector dark count rate determines the squeezing 
parameter in the virtual PDC's (see App. [Cj) . 
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FIG. 2: Given two copies of a single-rail entangled state p, 
dual-rail entanglement is obtained by mixing and conditioning 
on a single click at each end. The half-filled circles denote 
excitations shared between two modes. 



Memory imperfections (iii) are described by the phys- 
ical parameters of the particular quantum memory in- 
serted in the setup. However, it is desirable to analyse 
the problem using a general description of the memories. 
The incoming and outgoing light fields are conveniently 
described by harmonic oscillator degrees of freedom, and 
ideally the memories map the state of the incoming light 
to the outgoing light. In the Heisenberg picture this im- 
plies a mapping a' = a, where a' and a are the field 
operators for the outgoing and incoming modes respec- 
tively. In the presence of imperfections, the state transfer 
will be described by an admixture of other field operators 
into the outgoing mode. Here we are mainly interested 
in describing quantum memories based on atomic ensem- 
bles. Such ensembles can in the limit of many atoms be 
described by a set of harmonic oscillators. Furthermore, 
ensemble based quantum memories are to leading order 
typically described by interaction Hamiltonians which are 
quadratic in the field operators for the light and atoms 
0- E2; [HI- The resulting state evolution can then be 
described by a Bogoliubov (i.e. linear, unitary) trans- 
formation of the field operators even in the presence of 
imperfections such as spontaneous emission. We shall 
consider quantum memories described by the most gen- 
eral possible Bogoliubov transformation and our results 
thus apply to a very wide class of quantum memories. 
Our calculations cannot, however, describe interactions 
with Hamiltonians of higher order than quadratic in the 
mode operators. E.g. we cannot describe the optical 
Kerr effect or or single- atom memories. 

In addition to the errors (i),(ii),(iii) we also con- 
sider detectors which do not resolve the photon number, 
as efficient, single-photon counters are difficult to con- 
struct and it is interesting to compare the two cases of 
counting/non-counting. We note that the DLCZ setup is 
equivalent to the protocol considered here in the case of 
perfect memories, i.e. in the absence of type (iii) noise 
[2|. The results presented below therefore also apply to 
that protocol. 

The figures of merit for a repeater are the rate R at 
which entangled pairs are generated and the quality of 
the entanglement as functions of the distance L. As seen 
below, and as noted in Ref. for lossy or non-resolving 
detectors the final state will contain a large vacuum com- 
ponent with no excitations in the ensembles. However, by 
combining two single-rail qubits into one dual-rail qubit 
via postselection as shown in Fig. [2] the vacuum compo- 
nent can be removed and interesting applications such as 
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quantum cryptography or teleportation can be performed 
provided that the quality of entanglement conditioned on 
successful postselection is high [2, 

As a measure of entanglement, we use the Bell param- 
eter S. Let Pdiff(<p,tp) (P sa me(<f>, <fi)) denote the condi- 
tional probability that one upper and one lower detector 
(both upper or both lower detectors) on Fig.[2]click, given 
that a single click is detected at each end. We define 
E(4>, ip) = P S ame(<f>, ¥>) - Pdiff(<f>> <p) and 

IT IT 7T — 7T 7T — 7T 

S = E(-, -) + E(0, -) + E{Q, —) - E(-, —). (3) 

The parameter space is then — 2\/2 < S < 2\/ 2 with 
\S\ > 2 designating that the state is entangled [151 ]- Note 
that for ease of computation we have fixed the angles in 
the definition of S, whereas the custom definition is a 
maximum over varying angles. The angles in (J3j maxi- 
mize S for the ideal state p = and we prove 
(App. [Aj that to first order in perturbations away from 
the ideal state, fixing the angles gives the correct de- 
scription of the decrease in S. As an alternative to S one 
might calculate the fidelity of the postselected final state 
with respect to the ideal entangled state of two dual-rail 
qubits. The main goal of our analysis will be to deter- 
mine the dependence of S on the physical parameters in 
the setup. 

The rate is determined by the probabilities for suc- 
cessful entanglement generation, connection, and post- 
selection, denoted by qo, qi where i = l...n, and q ps 
respectively. The average time t„, it takes to create an 
entangled pair of length L, obeys the equation 



(4) 



where t' n is the average time it takes to create two neigh- 
bouring entangled pairs of length L and r = Lq/c is the 
classical communication time for an Lo-segment. The 
solution obtained by replacing t' n by t n in the above 
equation, i.e. by approximating the waiting time for two 
pairs by that of a single pair, considerably overestimates 
the rate. However, it turns out that a much better ap- 
proximation is found by moving this assumption one step 
lower, i.e. by instead taking the recurrence for t' n to be 



4+1 = ^+1(^2" + 4). 



(5) 



Here v n is the average number of tries needed for two in- 
dependent binomial events (i.e. generation/connection), 
each with probability q n , to both succeed 



_ 3 - 2q n 
(2-q n )q n 

and the recurrence has solution 

4 = r{2 n - x v n + ■■■ + 2% ...vx+v n ... vivo) 

M TV n . . . UlUQ, 



(6) 



where the last equality holds for qo <C n^ 1 since for any 
linear optical Bell measurement the success probability 
is at most 1/2 and therefore 
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(8) 



Since to obtain good entangled states it is necessary to 
keep the photon creation probability r 2 small, this condi- 
tion is well fulfilled in practice. Taking into account the 
final postselection step, we get the following expression 
for the rate 



R = t n l = t 1 q pa v n 1 



71+1 



q ps q n ■ ■ - qo, 



(9) 



(7) 



which is equivalent to the expression used in Refs. 0,0, 
0|. It also agrees with the empirical estimate found in 
Ref. Q. The last simplification is exact in the limit of 
small q and deviates from the first line of by at most 
a factor ~ 2.6 for n < 8. We have found that numerical 
simulation of the rate for given shows good agreement 
with 



III. METHODS 

To get at the figures of merit, we apply both analyt- 
ical and numerical methods. Our approach is to com- 
pute the two mode density matrix p n of the entangled 
pairs at each step of the protocol (i.e. as a function of 
L = 2 n Lo). Knowing the state, we may then calculate 
S and any other derived quantities. In the course of 
computing p n we also obtain the success probabilities qi 
for entanglement generation and connection, which de- 
termine the rate. Although the methods we have devel- 
oped are valid for arbitrary photon numbers, we shall in 
practise always work with small photon numbers, so that 
p n may be described by a 4x4 or 9x9 matrix in the Fock 
basis. 

Computing p n may be broken into two principal steps. 
First, computing the state po from entanglement gener- 
ation (Fig. Q] a) and second, computing p n for n > by 
iterating the connection process (Fig. [T]b). 

To deal with these two tasks, in particular the second 
one for which the Bogoliubov transformation depends 
on the atomic memory, we have developed a framework 
for calculating the output state from an arbitrary Bo- 
goliubov transformation followed by projective measure- 
ments, given the input state. Our method, described in 
App.[Bj is based on a generating function T . This func- 
tion takes two variables for each input and each output 
mode, and is defined such that its derivatives evaluated 
at zero form a matrix transforming the input to the out- 
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FIG. 3: Entanglement connection with no transmission losses. 
Two pairs of entangled light modes, each in state p n extending 
over a distance 2 n Lo, are connected. Detection of a single 
click heralds successful connection and leads to entanglement 
between the far-away ends of the pairs, leaving them in the 
state Pn+i extending a distance 2 n+1 Lo- A full state transfer 
through the memories is included in the connection step. 



put state, e.g., for a single input and output mode, 



\l\Pout 



\j) = ^Mijkl x (k\p in \l), 



(10) 



where unitarity requires that 



\h\ 



|ci| 2 -M 2 -|c 3 |^ = l. 



(13) 



Here the annihilation operator cti is for the input mode, 
and we have denoted the retrieved mode by a prime. The 
parameters bi,b 2 ,ci,c 2 , C3 are determined by the physical 
properties of the memory, and the state transfer is perfect 
when bi — 1 and only b\ is non-zero. Additional terms 
are due to memory imperfections. We note that with the 
proper choice of phases, all but one of the parameters in 
(fl2ll may be assumed real. See App.[C]for details. 

As an example of how we apply mode reduction and 
the generating function, consider now a very idealised 
repeater setup, where the errors are small and due solely 
to memory imperfections. We neglect transmission losses 
and we assume the entanglement generation to produce 
a perfect Bell-state, po = The parameters of 

the memories are for now assumed to fulfill 



b2,c 2 ,c 3 = 0, 



bf-1 



c\ < 1 



(14) 
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For any given Bogoliubov transformation and set of pro- 
jection operators, we can compute T and from T we find 
Pout for any given p in . 

In addition to the generating function, we also make 
use of mode reduction. The Bogoliubov transformation 
for a full state transfer (i.e. storage and subsequent re- 
trieval) of a light mode through a realistic atomic mem- 
ory often involves many auxiliary modes in addition to 
the input and output modes. However, in App. [C] we 
show that the number of modes can always be reduced 
to three. The most general transformation for the state 
transfer becomes 



a'i = b\a\ 



cia\ 



3202 



c 2 a3. 



C3S3' 



with b\, ci real. Using these assumptions we can com- 
( pute the generating function corresponding to entangle- 

ment connection and from the generating function we can 
find p n - It is convenient to include the full state trans- 
fer through the memories (light to atoms to light) in the 
connection step such that p n denotes the state of two 
entangled light modes after n connection steps. The con- 
nection then proceeds as shown in Fig. [3j The connecting 
ends of two neighbouring pairs of entangled light modes 
are stored in two memories, retrieved, mixed on a bal- 
anced beam splitter and then measured by non-counting 
photodetectors, with a single click heralding successful 
connection. The Bogoliubov transformation correspond- 
ing to connection may be found using lfl2|) and lfl4l) and 
T can then be computed from (|B4[) . Using T we find p n 
for the first few steps n = 1, 2, at each step expanding 
it to second order in c\. Based on the results we come 
(12) up with the following anzats: 
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(1 - 2/„)ci 
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- 9nC\ 
~ 9nC\ 




(1 - 2/ n )ci 



(1 - Vn)4 



(15) 



where p n is given in the Fock state basis, and f 7 i,g n are 
unknown, real valued functions to be determined. Our 
anzats will be confirmed, if the form (fl5]) is preserved 
under entanglement connection (Fig. [3|) and the result- 
ing recursion equations for /„, g n can be uniquely solved. 
Connecting two copies of (fl5| and expanding to second 



order in c\, we find that the form is indeed preserved if 

fn_ 

(16) 



fn+l — 



■Jn+1 



2 - fn ' 

4/n(4 + .g„) + H/„ 3 -20/„ 2 -4 



2/„(/„ - 2) 2 

For po to take the correct form, we must have fo = h, 
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go = 0, and these equations then yield the solution (see 
App.EJ 



fn 

9n = 



1 



2" + 1 ' 

_ 2 i+3« + 3 . 2 1+2n + 5 • 2" - 9 ' 17 ^ 
3-2(2™ + l) 2 ' 



Using L/Lq — 2" we arrive at the following result for the 
conditional Bell parameter: 



S = 2V2(1 - (L/L - lfc\) 



(18) 



to second order in c\. 

The above example illustrates the approach we take 
in deriving analytical results for the Bell parameter: Ex- 
cept for transmission losses in the detector arms, which 
may be considerable, errors are treated perturbatively 
and independently. First po is computed, and then an 
anzats for p n to the desired order in the error is found, 
leading to recurrence equations which are solved with 
initial conditions given by po- We treat the following 
errors perturbatively (in each case keeping transmission 
losses finite): finite initial squeezing, dark counts in en- 
tanglement generation, memory imperfections, and dark 
counts in entanglement connection. The results for the 
conditional Bell parameter in each case are presented in 
the next section. 

To verify our analytical results, we have also performed 
numerical simulations, where it is not necessary to treat 
the errors perturbatively or independently. We have com- 
puted the generating functions for the Bogoliubov trans- 
formations corresponding to Fig.[T](a), (b) and made use 
of these functions to numerically compute p n for vari- 
ous values of the losses, initial squeezing and dark count 
rates. Numerical results for repeaters using two specific 
atomic memories are presented in section Sec. IIV Bl and 
comparisons are made to the analytical approximations. 

The success probabilities qi in the rate are found by 
inserting the Bogoliubov transformations corresponding 
to entanglement generation or connection in (|B1|) and 
taking the trace. In this way one may in principle derive 
an expression for R valid for a repeater based on a general 
memory, as we have done for S. As we shall see below, 
however, the protocol only works well if the memories are 
close to being ideal apart from losses. We therefore only 
consider the rate for repeaters where losses are the sole 
errors since other imperfections will only slightly perturb 
the results. 



IV. RESULTS 

In the first section below we present our analytical re- 
sults for the conditional Bell parameter in the presence 
of errors, and in the following section the analytical ap- 
proximations are compared to numerical simulations of 
repeaters based on specific atomic quantum memories. 



A. Analytical results 

In the following subsections, we present results for 
S(L/L ) obtained by perturbation for each error source 
separately. Afterwards, we deal with cross terms between 
the perturbations. Arbitrary transmission losses are al- 
lowed in all cases and all results are given for both num- 
ber resolving and non-number resolving photodetectors. 

To get an intuitive idea of the nature of the errors, no- 
tice that the conditional Bell parameter S will decrease 
whenever an erroneous event during generation or con- 
nection, such as a dark count or memory imperfection, 
can lead to a non-vacuum, separable output. That is, 
errors occur whenever superfluous excitations are intro- 
duced into the system. Fig. shows how this may be 
caused by a connection dark count. Once an error has 
occurred, it propagates through the protocol as shown in 
Fig. HJd. Below we use such considerations to justify the 
scaling of S. 



Imperfect memories and connection dark counts 

Dark counts during entanglement connection can be 
treated as a memory imperfection by considering the vir- 
tual PDC and beam splitter introduced in Fig. Q] (b) to 
be a part of the memory protocol, and therefore we treat 
these two error sources simultaneously. An expression for 
the combined Bogoliubov transformation including dark 
counts is derived in App. [Cl 

We take the Bogoliubov transformation for a full state 
transfer through the atomic memories including dark 
counts to be (jClOP and neglect all other error sources 
except transmission losses. The photon loss probabil- 
ity in entanglement connection is taken to be p con - We 
treat each parameter in the transformation separately, 
and proceed as in the example given in Sec. IIIII One can 
check that to leading order, this perturbation contains no 
cross-terms, and hence we may add the errors in S deriv- 
ing from each parameter. For the case without photon 
counting, we find 



l-(L/L a -lj 2 (l+ Pcon ) 2 \ Cl \ 2 
-4(L/L ) 2 (l+p co „)(| C2 | 2 + |c 3 | i 



2^2 



(19) 



Tide 



1 -Pc 



and for the case with photon counting 



S 
2V2 



= l-(L/L Q ) 2 p 2 con \c 1 \ 2 (20) 



- 8(L/L ) 2 Pcon (\c 2 \ 2 + \c 3 \ 2 + ~r~^~~ ) 

^ Peon 



For the perturbation to be valid, the error in S must be 
less than 1, and so we require Ci, c%, C3 <C 1 and fide € 1- 
p CO n- In addition, here we consider the large L limit L Q <C 
L and give only the leading order in L. However, we have 
verified that the expressions are a good approximation 
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(a) 



|*+> 



(b) 




|*+) 



111) 



FIG. 4: Creation and propagation of errors. Filled circles de- 
note excitations, empty circles denote vacuum, (a) During 
connection of two segments in l^" 1 "), vacuum is read out but 
the connection is accepted due to a detector dark count leav- 
ing the remaining modes in the separable state |11). (b) Once 
an error has occurred, it may propagate. Connecting 1 11) with 
and ideal entangled state and requiring a single click leads to 
1 11) in the output (if no loss occurs). A vacuum state of the 
form 1 00) propagates in the same fashion. 



to the exact analytical result for S (which could not be 
put on a closed form) also for L ~ L , as long as the 
perturbative condition is fulfilled. 

There are several things to notice about the results lfT9|) 
and (|20p. First, note that the conditional Bell parameter 
is independent of the parameter 62. This is because 62 
corresponds to a plain loss. If 62 is the only imperfection, 
the transformation fT2|) is passive and hence 62 leads to 
an increase of the vacuum component of p n , which does 
not influence the conditional Bell parameter. For the 
same reason, there is no term in S depending only on 
the transmission loss p CO n- Second, note that the errors 
are suppressed for vanishing p CO n, when the photons are 
counted, but persist for vanishing p con , if they are not 
counted. This, and also the scaling of the error with 
L and p con , may be motivated by the following simple 
picture. 

In the case of only c\ non-zero the unitarity condition 
(fl~3j) becomes |&i| 2 — \c\\ 2 = 1, and it follows that the 
effect of non-zero c\ is a one-mode squeezing of the in- 
put mode. Similarly the effect of non-zero C2 or C3 is a 
two-mode squeezing of the input mode with an auxiliary 
mode. Letting p = p n <£> p n the memories take 



or 



P 



P 



• Si(d)S fl (ci)pS£(ci)£k(ci) (21) 



SLl{c2)SR r (c 2 )pSt,(c2)S' 



Rr 



(<*) 



(22) 



where (referring to Fig. [3j L, R are the measured modes, 
I, r are two auxiliary modes and Si(c) = exp(ca| 2 — c*af) 
and Sij(c) = exp(caja] — c*ajdj) are the single and two- 
mode squeezing operators. From these expressions we 
see, that ci errors introduce photon pairs into the mea- 
sured modes, while c 2 and c 3 errors introduce single pho- 
tons. To lowest order in the c's the amplitudes for errors 
to occur are c\ , C2, C3 and hence the error in S scales 



with I ci 1 2 , 1 C2 1 2 , 1 C3 1 2 . Now, the post-selection implied by 
S (see Fig. [2]) requires that any superfluous photons be 
removed before the post-selection stage. In the case of 
photon counting detectors, this can happen only through 
photon loss, and therefore the ci and 02,03 error terms 
scale with p 2 con and p con respectively. In the case of non- 
resolving detectors, in addition to loss, superfluous pho- 
tons can be removed when multiple photons are incident 
on the same detector producing only a single click. This 
implies that errors persists for vanishing loss, and is ap- 
parent by the replacement p con — » (1 +p con )/2 from lfl9|) 
to <f20|) . To understand the scaling with L, note that there 
are L/Lq — 1 connection attempts in total. There are 
therefore L/Lo ways for a photon to get lost and L/Lq 
ways for an extra photon to be introduced. In a C2 or C3 
error one extra photon is introduced and thus one pho- 
ton must get lost, whereas in a c\ error two photons are 
added and two photon must get lost. However one of the 
added photons must get lost in the connection attempt in 
which it was created, since successful connection requires 
exactly one detector click, and the scaling for both types 
of error is therefore L 2 . Since dark counts can be treated 
by mixing of the signal mode with a two mode squeezed 
state (see Fig. [1]) this error term scales the same way as 
the C3 term. This is also apparent from the Bogoliubov 
transformation (|C10|1 . Note that as p CO n — > 1 , the proba- 
bility that a given click is a dark count approaches 1 and 
the dark-count error term diverges. 



Finite initial squeezing 

Ideally, each entangled state consists of one excitation 
shared between two modes, and as seen above the in- 
troduction of additional excitations is a cause for errors. 
For this reason, when non-counting detectors are used 
or loss is present, the squeezing parameter of the state 
(HJ) used for entanglement generation must be small and 
ideal Bell-type entanglement is only achieved in the limit 
of negligible squeezing r — > 0. Of course this ideal limit 
for the Bell parameter is the worst-case limit for the rate, 
as the success probability for entanglement generation is 
proportional to the probability r 2 of creating a photon 
pair in the PDC's. Thus, there is a trade-off between 
the conditional Bell parameter and the rate, and we are 
required to keep r finite. Here we examine the effect of 
finite initial squeezing on S for the final entangled state 

Pn- 

The state po produced by the setup Fig. Q] (a) can be 
calculated either directly in the language of mode opera- 
tors starting from the state (Q} or by computing the gen- 
erating function of the entire setup including the PDC's 
and assuming vacuum input. We have found p n to second 
order in r, neglecting dark counts but making no restric- 
tion on the loss parameters p ge n,Pcon- Going to second 
order in r is equivalent to taking a maximal photon num- 
ber of two, and hence po is described by a 9x9 matrix in 
this case. Knowing po we proceed as in Sec. lIIII to find p n , 
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assuming perfect entanglement connection. The results 
for the conditional Bell parameter are 



S 
2V5 



1 - 8(L/L ) 2 1 + o Pge " 1 + f °" r 2 (23) 



without photon counting, and 



2\/2 



= 1 - 8(L/L ) PgenPcon T 



(24) 



with photon counting. 

The factor of r 2 indicates that, similar to what we saw 
previously, the decay of S is caused by superfluous excita- 
tions in the repeater, in this case coming from the PDC's. 
The extra photons come in pairs, with one in the detec- 
tor and one in the memory arm of Fig. Q] (a), and hence 
the situation is analogous to that of ci-type errors above. 
Two photons must get lost before the post-selection stage, 
and one of these must get lost in the generation attempt 
in which extra photons were introduced, since a single 
detector click is required for successful generation. With 
photon counting this type of error is suppressed when 
there is no loss in generation (p gen — > 0) or connection 
(Peon — + 0) , but persist for vanishing loss if the photons 
are not counted, and the error terms scales quadratically 
in L. 



Generation dark counts 

Last, we consider the effect of dark counts in the entan- 
glement generation, but neglect errors of second order or 
higher in r. Again, the generated state can be derived by 
means of the generating function for the setup Fig.[T](a). 
Having derived po we proceed to find p n and the con- 
ditional Bell parameter assuming perfect entanglement 
connection. To lowest order in the dark count probabil- 
ity fide we find, without photon counting 



-^= = 1-4(L/L ) 2 i±te 

2V2 1 - Pgen 



ndc, 



and with counting 
S 



2V2 



1 - 8(L/L ) 2 



Pge 



1 



Pge 



-ndc 



(25) 



(26) 



Again, notice that the error in S is suppressed for 
Pgen — » in the counting case but not in the non-counting 
case. Because a dark count alone results in the generation 
of a vacuum state, it must be combined with a double ex- 
citation in some segment to lead to an error in S which 
explains the quadratic scaling of the error. Despite the 
fact that a double excitation is needed, r 2 does not show 
up in eqns. (|25|) . l(26|) . because although the probability 
for a double excitation is of order r 4 , there is no need to 
create an excitation in the segment where the dark count 
occurs, and hence the total number of generated photons 
is unchanged with respect to the ideal case. 



Perturbation cross terms 

So far we have considered each error source indepen- 
dently, and we have found perturbatively their effect on 
the conditional Bell parameter. However, we have not 
addressed the possibility for cross terms in the perturba- 
tion, when errors of different type are present simultane- 
ously, as will always be the case experimentally. Indeed, 
cross terms do appear and may have a severe effect on 
S for certain values of the parameters. Here we identify 
the regime where cross terms can be safely neglected. 

We find that the significant cross terms are those aris- 
ing from the combination of a generation dark count with 
a memory imperfection or a connection dark count. To 
see this note, that a generation dark count results in 
a vacuum state, and connecting this with a Bell state 
in the presence of a memory imperfection or a gen- 
eration dark count may result in the separable state 
1 01) , which leads to an error in S. Since the event re- 
quires a generation dark count and an error in connec- 
tion, the error term must be proportional to ehd c , with 
e = | ci | 2 , | C2 1 2 , | C3 1 2 , or n^c- However, compared with 
the errors considered in previous sections we now require 
generation of only L/Lo — 1 photons rather than L/Lq, 
and therefore this error term is also enhanced by a fac- 
tor of 1/r 2 , so that the total order of magnitude of the 
term is efidc/r 2 . Relative to the error terms considered 
previously, there is a factor of fidc/r 2 which amounts to 
enhancement if the dark count probability is higher than 
the probability to generate a photon, and to suppression 
in the reverse case. On the other hand, the cross term in- 
volving a generation dark count and an additional gener- 
ated photon, or the cross terms not involving generation 
dark counts must have L/Lq generated photons. They 
are therefore of order e 2 and can be safely neglected in 
both cases. 

We conclude that the analytical results derived in the 
previous sections provide a full description of the condi- 
tional Bell parameter whenever the production rate for 
photon pairs in the parametric down conversion is signif- 
icantly higher than the dark count rate of the detectors. 
Since the errors due to finite initial squeezing and dark 
counts scale in roughly the same manner according to 
Eqns. Hal, PD]1. p3 > p> ]) . it will be advantageous to 
make r 2 at least comparable with fide to increase the 
rate. We are then only making a minor error by neglect- 
ing cross terms. 



B. Application to specific memories 

We now verify the analytical results of the previous 
section by comparing with numerical simulations. The 
results for dark count and initial squeezing errors are in- 
dependent of the memories and hence apply to any DLCZ 
repeater with the architecture presented above. We re- 
mind the reader that this includes the original DLCZ 
protocol since the interaction used by DLCZ for entangle- 



2" 2 2 o 2 2 o 2 2 

Distance L/Lq 



2" 2 2 2 2 2 2 2 2 

Distance L/Lq 



FIG. 5: The effect of finite initial squeezing. S(L) is plotted 
for r = 10~ 2 and transmission losses of p a en — U.<y, peon — 
0.1, in the case of counting (dots) and non-counting (circles) 
detectors. Dashed lines show the analytical approximations. 



ment generation is effectively a two-mode squeezing and 
hence equivalent to our PDC's 0|. For the memory re- 
sults we perform simulations using two specific ensemble 
based quantum memories. One proposed by Muschik and 
Hammerer [12], which we denote the two-pass memory, 
and one proposed by Julsgaard et al. here denoted 
the one-pass memory. 

As explained in Sec. Hill our numerical results are ob- 
tained by computing p n by means of the generating func- 
tion for specific values of the parameters of the system. 
In doing so, there is no need for treating errors pertur- 
batively but we do have to restrict the dimension of p n , 
because it is not practical to work with large matrices. 
Since p n is given in the Fock state basis, this effectively 
means restricting the maximal number of excitations in 
the entangled states. In our simulations we take the max- 
imal excitation number to be 2 so that p n is 9x9. This 
implies that our numerical results can be considered ex- 
act as long as the probability of creating three or more 
excitations is negligible, however for the repeater to work 
reliably beyond a few connection steps r 2 <C 1 is required, 
so this condition is well fulfilled in practice. 

Results for the effects of finite initial squeezing and of 
dark counts on the conditional Bell parameter are dis- 
played in Figs. I5|6l For each type of error we display 
the analytical conditional Bell parameter from the pre- 
vious section as well as the numerical simulation. Note 
that in all cases, S is well described by the analytical 
approximations for distances L where S is well above the 
classical threshold of 2. The excellent agreement con- 
firms the analytical results of Sec. IIV Al The perturba- 
tive approximation breaks down at the distance where S 
drop below the classical threshold of 2, i.e. the maximal 
distance over which communication is possible with the 
given losses, since at this point the error in S is no longer 
small. 

We can estimate the maximal distance imposed by 
dark counts or PDC squeezing by setting S — 2 in f23|) 



FIG. 6: Effect of detector dark counts in entanglement con- 
nection. S(L) is plotted for fide = 10~ 4 and transmission 
losses of Pg e n = 0.9, peon = 0.1, for counting (dots) and 
non-counting (circles) detectors. Dashed lines are analytical 
results. The inset shows the effect of dark counts in entangle- 
ment generation, for the same parameter values. 



and |(25]) to get 



(L/L ) 



V2 



(L/Lo) 2 



4(1 +Pgen)(l + Peon ) 

2- V21 



8 



1 



Pge 



qen 1 



(27) 
(28) 



in the non-counting case. Small r are not hard to ob- 
tain experimentally, but the dark count number fide is 
dictated by the detector dark count rate and the pulse 
duration r. With microsecond pulses, an optimistic 1Hz 
dark count rate, and 90% photon loss in generation, the 
maximal distance ([28| is L ~ 64L . Recent protocols 
offer better tolerance for dark counts and other multi- 
photon errors, see e.g. Ref. 

As mentioned in the previous section, the different er- 
rors introducing superfluous excitations in the repeater 
can only be treated independently in perturbation theory 
when the production rate in the PDC's is higher than 
the rate of spurious excitations. This is demonstrated 
in Fig. [3 where we plot S in the two cases fide < r 2 
and fide > r when both generation and connection dark 
counts occur. We see that the analytical approximation 
obtained by adding the error terms in lfl9]) and (|25| is 
valid only in the former case 

Next we turn to specific quantum memories. The two- 
pass memory is depicted in Fig. [8l It consists of an en- 
semble of atoms contained in a glass cell at room tem- 
perature and traversed twice in orthogonal directions by 
a light beam. Both light and atomic spin are strongly 
polarised, such that the Holstein-Primakoff approxima- 
tion may be applied and each system can be described 
by a harmonic oscillator. The interaction Hamiltonian 
is xaXl in one pass of the light pulse and paPl in the 
other, where subscripts A and L refer to atoms and 
light, and the overall interaction Hamiltonian becomes 
XaXl + VaVl ~ a\a-L + o>ao>l- As a resu lt of this inter- 
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FIG. 7: Significance of perturbation cross terms. S(L) is plot- 
ted with dark counts in all detectors for fia c = 10~ 5 in the two 
cases r 2 = 10" 4 (left )and r 2 = 10~ 6 (right). When n dc > r 2 , 
a large discrepancy is observed between the numerical result 
and the analytical approximation omitting perturbation cross 
terms (dashed line). 



action, the polarisation state of the light and the state of 
the collective atomic spin are swapped, and this process is 
governed by the light-atom interaction strength k. In the 
absence of losses a full state transfer through the memory 
may be described by the Bogoliubov transformation (see 
Ref. ; 1 2| for details) 



(e-* 2 - l)a L - e-« 2 / 2 VT 



e- K / 2 aT- 



(29) 

where a,L,a' L are the stored and retrieved light modes, &a 
is a collective atomic mode and a r £ l is the input mode 
of the retrieval light pulse. From i|29p it is apparent that 
the memory is perfect when k is large. However, in ex- 
periment k is restricted by spontaneous emission, and 
in the following examples we use an optimistic value of 
k = 2. In the lossless case |29|) is a passive transforma- 
tion and hence does not introduce any error in S (c.f. 
(H2j) , (Unj) , (EOj ) . Only the rate is affected by k. 

The important error parameter in the setup Fig. [8] is 
the reflection coefficient £ for reflections at the cell walls. 
Reflections occurring between the two passes of the light 
pulse introduce an active part to the transformation f29|) 
leading to non-zero C3 in (fT2| . and in fact we find cf « 
0.9£ to lowest order in £ when k — 2 [16]. Putting 5 = 2 
in lfl9l) we have 



V2 



0.9-8(l+p co „) 



(L/L ) 



(30) 



Taking an optimistic value £ = 10~ 3 and no loss, this 
limits the communication distance to L ~ 8Lo- In Fig. [9] 
we plot S for the two-pass memory together with our 
numerical approximation. 

The one-pass memory has been demonstrated experi- 
mentally by Julsgaard et al. [13]. The write-in setup is 
similar to the two-pass setup in Fig. [8] but instead of the 
polariser a detector is placed is placed in the beam and 
the light is measured after one pass through the atomic 
ensemble. A feedback is then supplied to the atoms via 
a magnetic pulse based on the measurement outcome. In 




FIG. 8: Setup for the two-pass memory. An atomic ensemble 
is placed in a magnetic field, with the spin of the ensemble 
strongly polarised along the field. The atoms are traversed 
twice in orthogonal directions in the plane perpendicular to 
the field by a light beam. As a result, the light polarisation 
state can be stored in or retrieved from the atomic spin. 
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FIG. 9: Effect of two-pass memory imperfections. S(L) is 
plotted for £ = 10~ 4 , 10~ 3 and transmission losses of p ge n = 
0.9, Peon = 0.1, for counting (dots) and non-counting (circles) 
detectors. Dashed lines are analytical results. 



order to see whether this memory is suitable for use in 
the repeater, we simply assume perfect readout. With 
this assumption the Bogoliubov transformation of a full 
state transfer becomes 



Kg . „ Kg „ + 1. . „ 1. .+ 
a L = (1- -Y)aA + -^a> A + -(K+g)a L --{K-g)a^ (31) 



where k is the coupling strength and g is the feedback 
gain. From this expression it is clear that the memory is 
never perfect, since for any non-zero choice of k and g the 
output light contains some mixing in of the atomic mode. 
However, if the atomic mode is squeezed prior to storage, 
the noise introduced by the atoms can be suppressed. 
Assuming that the variance in the A"-quadrature of the 
atomic mode is squeezed by a factor s, and putting k = 



10 



3n 




Distance L/Lq 



FIG. 10: Effect of one-pass memory imperfections. S(L) 
is plotted for squeezing parameters (from the top down) 
s = 10 5 , 10 3 , 10 and transmission losses of p gen = 0.9, p CO n = 
0.1 for counting (dots) and non-counting (circles) detectors. 
Dashed lines are analytical results. 

3=1 the Bogoliubov of a state transfer becomes 

a' L = ^L(a A + a\)+a L . (32) 

This is on the form lfT2]l and it is now easy to read of the 
coefficients and plug into lfl9|) . Taking S — 2 and solving 
for s we get 

s=^ ^^- (V^o)- 2 ~-30dB (33) 

£ J- "T Peon 

for L = 16Lo and no loss. This value of s is far beyond 
what can be achieved experimentally at the moment [13l . 
[13, Hi (3- E.g. in one recent experiment about 3dB of 
squeezing was reported fl9| . We therefore conclude that 
the one-pass memory is not suitable for implementation 
of a DLCZ-type repeater. The performance for several 
values of the squeezing is shown in Fig. [10] 

Finally, let us briefly consider the rate. We prove in 
App. 10 that the rate of the DLCZ-repeater in the ab- 
sence of dark counts and memory imperfections is given 

by 

R= ^r 2 {l-p gm ){L/L Q )-^l 2 R> (34) 

with 

R > = e jhMi^^ L / L °)-ThM P T^^). (35) 
in the non-counting case, and 

R=^r 2 (l-p gen )(L/L r lo ^ 3 R / (36) 

with 

R > = e TTr2Li2( 5 ^ T V2L/Lo)- T ^Li 2 ( ^^jV?) (37) 

in the counting case. We can obtain the scaling of the 
rate for a fixed final imperfection in S , by using 1)231) or 
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FIG. 11: The rate for a fixed imperfection in S of 5%. We plot 
the two cases of counting (circles) and non-counting (dots) 
detectors, for p gen = 0.9 and p CO n = 0.1 (upper curves), 
Peon = 0.9 (lower curves). For reference we plot the ideal 
scaling L~ 2 - 10 ^ 3 (dashed line). 



([24"ll to determine the value of r and inserting this value 
into the above expressions. Fig. Qj] shows the result. In 
the best case, when the connection loss is negligible, the 
scaling is L _2_lo S2 3 _ When connection losses are small, 
the rate is seen to be significantly enhanced by the use of 
counting detectors. This can be understood as a conse- 
quence of the vacuum component of the state p n growing 
faster with n for non-counting detectors. Connecting two 
entangled pairs using non-counting detectors may lead to 
a vacuum state even in the absence of any losses. This 
is not the case for counting detectors. A fast growing 
vacuum component leads to a low connection probability, 
since the vacuum cannot contribute to the clicks required 
for successful connection. As losses increase, the proba- 
bility for two photons to reach a detector simultaneously 
decreases, and hence the advantage of counting over non- 
counting detectors disappears, the two rates being equal 
in the limit of very high loss. 



V. CONCLUSION 

We have presented a thorough analysis of quantum re- 
peaters using general ensemble based memories in the 
DLCZ architecture. As a primary result of our analysis 
we have derived perturbative analytical expressions for 
the Bell parameter of the generated entangled states in 
terms of the distance and the memory parameters. Our 
results apply to repeaters based on any quantum memory 
which may be described by a Bogoliubov transformation 
and hence to any system for which the interaction is no 
more than quadratic in the mode operators. We have 
verified our analytical results by comparison to numeri- 
cal simulations and found good agreement in the range 
where perturbation theory applies. 
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We are aware that many protocols improving upon 
the DLCZ repeater architecture have been put forward, 
which promise significantly better tolerance for multi- 
photon errors and better rates, and that hence the system 
analysed in this paper is not a likely candidate for experi- 
mental implementation @, 0, S, S E3| ■ However, we have 
demonstrated how to analyse a repeater in terms of a very 
general quantum memory model, and our methods can 
in principle be extended to any system for which entan- 
glement generation and creation are described by Bogoli- 
ubov transformations and single-photon detection. This 
is the case for the recent proposals Refs. [!, 0, [H, [f| EH • 



where A 



6? and we have used that 4r- 



0,0° 



0. 



Thus we have proven that to first order in the perturba- 
tion, it is optimal to evaluate S using the angles cp . 



APPENDIX B: THE GENERATING FUNCTION 
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APPENDIX A: BELL PARAMETER ANGLES 



We define the generating function in the general case, 
where our system S has arbitrarily many modes, some of 
which are output while the remaining modes are mea- 
sured or traced out. Symbolically S = OIZ, where O 
are the output and 1Z the remaining modes. If the sys- 
tem is subject to a Bogoliubov transformation U and a 
subsequent measurement with an outcome corresponding 
to projection operator P, then the unnormalised output 
state is 



Here we prove that the angles in the definition of S can 
be held fixed when calculating the decrease in S in first 
order perturbation theory. Let cp° = (0°, . . . , 0°) be the 
angles maximising 5(1 ^P" 4 ") ( V I /+ |) and let p(x) be a pertur- 
bation away from the Bell state with p(0) — |^ r+ )(^ r+ |. 
Then the angles which maximise S(p) will also be close 
to <p a and to first order in the perturbation, we can write 



S(x, cp) = 2V2 + J2 



i=l 
k 



»=i 

= 5(x,0°), 



dS 

dxi 

dS 

dxi 



o,4>" 



0,0° 



E 



dS 



0,0° 



(Al) 



P°o ut = TTn(PUpftfpi). (Bl) 
To construct the generating function, first note that a 
Fock state in mode i may be written 



= -^=(a,t)"|0) i = 



1 d n -t 



(B2) 



a=0 



where the parameter a can be chosen real and |0)i de- 
notes the vacuum state of mode i. The generating func- 
tion is defined to be 



•F(ai,/3i,7i,5i) = Tr K 



e 2 x (vac 



o(0| [j[e s >*°jPU \J[M 



vac)(vac\ ^J^e^^J C/ t P t ^Q e7o<i °j l°>o 

u ] p ] (j{D { lo )\Q) {()\D {-5 )^ pu njA(A) 
I 



(B3) 

vac) (B4) 



where i runs over all of S and o runs over the output sion (|B4|) with S = J2i( a i + Pi) + So(7o + $0) follows 
modes O. The parameters <Xi,/3i and ji,Si are real and from the circular property of the trace and the relation 
correspond to the input and output modes respectively. 

Using (EH) and (B3} we can generalise JTSJ, and we see e Q&+ 10> = e^l a l 2 D(a)|0). (B5) 

that T generates the output state as desired. The expres- 

which is a consequence of the disentangling theorem [20(]. 
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We now show how to obtain the function T for given 
U and P. Let the Bogoliubov transformation U be given 

by 

tftfij 17 = b a & i + c a 4 ■ ( B6 ) 

i 

To compute T we make use of several properties of dis- 
placement operators. First, one may prove that the vac- 
uum projection operator can be written as the integral 

|0)(0|=/^e-(^)/^(^). (B7) 

It follows from l|B2j) that the projection on any Fock state 
can be written in terms of derivatives of an integral over 
a product of displacement operators. Second, under U 
the displacement operators transform as 

rfDjifiU = J] A(/3^ - (3* Cji ). (B8) 

i 

And third, the product and vacuum expectation value of 
displacement operators are given by 

D{a)D{f3) = e* Im ( Q/r )£>(a + /?), (B9) 

(0|£>(a)|0) - e^l"l 2 . (BIO) 

Starting from (|B4|1 . the function T is found in four 
steps. First, all projection operators in the expression 
are replaced by integrals of displacement operators by 
making use of (|B7|1 . Second, the Bogoliubov transforma- 
tion is eliminated from the expression via l|B8p . Third, 
the integrals are pulled outside the vacuum expectation 
which then contains only a product of displacement op- 
erators, and the expectation value is evaluated by using 
(|B9|I and (|B10|) , Last, the resulting Gaussian function 
is integrated and we obtain an analytic expression for T 
involving only the a, (3, 7, (5-variables and the parameters 
of U. 



Projection operators 

The generating function as defined above can be com- 
puted for any measurement described by a projection in 
the Fock state basis. However, in this article we are in- 
terested only in measurements where a single or no click 
is observed. The measurement operator corresponding to 
the zero outcome (no click) is simply a projection on the 
vacuum state of the measured mode: 

7 J darfe = |0)(0|. (Bll) 

The measurement operator corresponding to a single 
click depends on the resolution properties of the detec- 
tor. We work with two contrasting cases. For perfect 
single-photon counters, the operator is 

*W,t = |l>(l| (B12) 

= ^L e (l Q l 2+ l' 9 l 2 )/ 2 /J(a)|0)<0|/J(/3)l 



where 1|B2|1 and <|B5|) was used. For detectors which can 
only distinguish between presence and absence of light 
but give no information about the photon number, the 
operator becomes 

Plight = 1- |0)(0|. (B13) 

Squeezed initial states 

In principle the generating function l|B3P allows for any 
input and output states. In practice, it is not conve- 
nient to work with high photon numbers, because the 
density matrices become large and the calculation of the 
elements corresponding to many-photon Fock states re- 
quire derivatives of high orders (c.f. (JTUJ) ) - One might 
therefore expect that squeezed initial states would be 
treated only approximately, as they have non-zero over- 
lap with all the Fock states. Fortunately, the squeezed 
states belong to a class of initial states which may be 
treated exactly: For any state, which can be written on 
the form A\vac), where the operator A generates a Bo- 
goliubov transformation, we see from <|B1[) that we may 
replace U by UA and take p l g = \vac)(vac\. In partic- 
ular, for the two-mode squeezed state ([I]) in modes 
we let U — > USij (r) where S# (r) is the usual two-mode 
squeezing operator. In this way the input is treated ex- 
actly, regardless of photon number. 

APPENDIX C: MODE REDUCTION 

In this appendix we show how to reduce the most gen- 
eral Bogoliubov transformation to three modes, and how 
to include dark counts in the memory Bogoliubov trans- 
formation. The most general Bogoliubov transformation 
for a full state transfer through a quantum memory takes 
the form 

= bi&i + c\a\ + bjCLj + Cja\ (CI) 

i 

where Si is the mode operator of the stored mode, and 
a-i that of the retrieved mode. Now, it is always possible 
to define new, independent mode operators 02,0.3 by 

b 2 a 2 = ^2 Mi, c 2 a\ + c 3 a| = ^ CjoJ, (C2) 

i i 

where &2,C2,C3 are complex coefficients, and hence we 
can simplify (|C 1[) : 

a'i = biai + c\a\ + b 2 a,2 + c%a\ + c 3 a 3 . (C3) 

We see that the new transformation involves only three 
modes. 

The b,c coefficients are determined by (|C2[) and the 
canonical commutator relations for a 2 ,a 3 . There is some 
freedom in the choice of phases however. The phases of 
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bi, ci can be adjusted by simple phase shifts of the input 
and output modes. It can be seen, that in our repeater 
setup a phase shift on the output mode has no effect 
on the measurement outcomes when connecting pairs, 
and hence we may always assume that either b\ or c\ 
is real. It is not obvious that choosing both of them real 
corresponds to an optimal phase choice in terms of S, 
however we have checked numerically that a phase change 
of the input mode has negligible effect on S. Hence in 
this paper we take b\ and C\ to be real. Any complex 
phase on b 2 can be absorbed into the definition of d 2 , 
and likewise the phase of C3 can be absorbed in 0,3. We 
may therefore assume b 2 , C3 to be real. Putting things 
together, we have 

h,b 2 ,ci,c 3 €R, c 2 eC. (C4) 

And because ai must preserve the canonical commuta- 
tion relations 

b\ + bl-cl-\c 2 \ 2 -4 = 1. (C5) 

Including dark counts 

As an example of mode reduction, we consider dark 
counts in entanglement connection. They are treated by 
including a PDC and beam splitter after the memory 
readout, as shown in Fig.[T](b). Tracing out one mode of 
the PDC leaves a thermal state in the other mode, which 
is mixed with the signal from the memory. We now show 
how the memory transformation l|C3P can be modified to 
take account of the dark counts. 

The new Bogoliubov transformation, including the ad- 
ditional modes a sl ,a s2 , is given by U B sU m e m S s i S 2{r), 
where Ubs is a beam splitter transformation, U mem is the 
transformation <|C3j) and we have included the squeezing 
as explained in App.[Bj This implies 

o'i = Sl ls2 (r)(af 1 -y/T^p + a sl ^/p)S sls2 (r) 

— a\y/\ — p + y/p (fi s i coshs — a\ 2 sinh s). (C6) 

Using l|C3j) and applying mode reduction, we can write 

a" = b'xai + c'^l + b' 2 a! 2 + c' 2 a$ + c' 3 a] , (C7) 

where 

b 2 a 2 = y/l -pb 2 a 2 + ^/pcosh(s) d s i, (C8) 
c 2«2 + c 3«3 = V 1 -P{ c 2a 2 + c 3 &l) - ^/psinh(s)al 2 . 

Choosing 61',62',c3' real and positive and using the 
canonical commutators, we find the primed coefficients 
to be 



b[ = 


yfl -pbi 


b' 2 = 


J (1 — v)b\ + pcosh(s) 2 


c[ = 


V 1 -pci 


4 = 


(1 - P)b 2 c 2 /b' 2 


4 = 


J(l-p)(\c 2 \2 + 4)+p S inh(s)*-\c> 2 \2 



It remains to relate the parameters p, s of the virtual 
optical elements to the physical dark count rate. One 
may prove, that the output state from one arm of the 
PDC, when the other arm is traced out, is a thermal state 
of mean photon number sinh(s) 2 . The average number 
of dark counts must equal the mean photon number at 
the detector due to the virtual PDC, and we therefore 
get fide = psinh(s) 2 . To obtain a final expression for 
the Bogoliubov transformation including dark counts we 
rewrite <|C9[1 in terms of the physical parameter fide, and 
since we are only interested in introducing dark counts 
but not photon loss to the memory output mode we let 
p — * while keeping fide constant. The result is 

b[ = h 

b' 2 = \J b\ + fide 

ci = ci (CIO) 
c' 2 = b 2 c 2 /b 2 

c' 3 = \J\C 2 \ 2 - \c' 2 \ 2 + C 2 + fide- 

Note that since we have chosen to include the dark 
counts in the memory transformation, attention should 
be payed to keeping the dark count rate fixed when pho- 
ton loss (peon > 0) is introduced. 

APPENDIX D: MATHEMATICAL DETAILS OF 
THE RECURRENCE EXAMPLE 

To solve the recurrence equations ([I6)l . make the vari- 
able substitution /„ = Z" 1 to obtain: 

fn+l = 2/„ - 1. (Dl) 

This equation is easily solved, subject to the initial con- 
dition fo = f^ 1 = 2, and we find: 

In = 2" + 1. (D2) 

Inserting the solution for /„ into the g„-recurrence and 
making the substitution g n = 2(2™ + l)g n one finds the 
recurrence equation 

g n+1 = 2g n - 2 3n+2 + 2 2n+2 - 3. (D3) 

Given the initial condition g = 2(2™ + l)g = 0, we get 
the solution 

g„ = —(2 n - l)(2 2n+1 - 2"+ 2 - 9). (D4) 
3 

From (|D2|) . l|D4|) and the definitions of f n , g n we have the 
solutions lfl7|) . 

Having found /„, g n we may obtain the conditional Bell 
parameter of the state l|15p as a function of entanglement 
distance. The conditional Bell parameter of the state is 

S = 2V2(f n - c 2 g n ) 2 /[fl - (2f n9n - (2/„ - l) 2 )c 2 

-((2/„-l + .g n ) 2 -2 5 2 )cf]. (D5) 
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FIG. 12: Plot of the rel ative error (LHS - RHS)/LHS where 
LHS and RHS refer to pTT)t . 



Now, if ci is small such that /„ dominates the denomi- 
nator, we can plug in the solutions for /„, g n and expand 
to lowest order in c\ . This gives the expression (fl8| . 

The derivation of the reduction in Bell parameter in 
the presence of fibre loss, and due to other error sources 
(other memory imperfections, dark counts and finite ini- 
tial squeezing) proceed along the same lines as the deriva- 
tion presented here and in Sec. IIIII However, the re- 
currence equations tend to be considerably more com- 
plicated when losses are included, and in some cases we 
have not been able to obtain a closed form analytical so- 
lution. In those cases we have obtained an exact solution 
of the recurrence numerically (by substituting the equa- 
tion into itself) and from this solution we have deduced 
the behaviour at large L/Lq. Subsequently we have ver- 
ified, by comparing with numerical simulations, that the 
analytical expressions thus obtained are also valid for L 
close to La. 



APPENDIX E: THE DLCZ-REPEATER RATE 
WITH ARBITRARY LOSSES 

We consider entanglement generation and connection 
as shown in Fig. [H and concentrate on the case where 
the photon production rate in generation is very low 
and where the memories are passive, such that no multi- 
photon errors are present. For such a system, entan- 
glement generation produces the state po = \^ + }(^ + \. 
Since photon loss is the only error, we expect p n to take 
the form 



Pn = Vn\^ + ){^ + \ + (1 - r} n )\vac)(vac\, 



(El) 



where i] n is a number and 770 = 1. If this form of p n 
is conserved under entanglement connection, Fig. Q] (b), 
then it follows by induction that it is correct. Assuming 
that the detectors resolve single photons, it is not difficult 



to see that this is indeed the case, and that 

tin 



Vn+1 = 

which has the solution 

Vn = 



1 



(E2) 



(E3) 



1 Peon ~t~ ^ n Pcon 

Now, from IjEip the success probability for connection is 



Qn+l = ^(1 -Pcon)l]n(2 - 7? n (l ~ Pcon)) 



1 



(1 -Pcon)vl/Vn+l 



(E4) 



Since r <C 1, the success probability for generation is qo = 
2r 2 (\—p gen ). The probability for successful postselection 
is q ps = Vn/2 an d from |(9|) the rate is then 

2 " 
R = ^TT 7 ^ 1 -JWX 1 ~ PconTvt^- (E5) 



We now put this expression on a closed form. We start by 
turning the product into a sum by taking the logarithm 

n n 

in n w = ~ n in ( 1 ~ pc ° n + Tpc ° n ) • ( Ee ) 

i=l i—1 

The sum can be estimated by taking the integral 

r n+l 

In (1 - Pcon + p conl 2 x )dx, (E7) 



where we have introduced a constant 7. By adjusting 7 
we make sure that the integral agrees with the sum above 
in the limits where the sum can be easily evaluated. The 
integral gives 



nln(l -p co 



Li: 



In 2 



(E8) 



where IA2 is the dilogarithm [21(. Since Lia(O) is 0, the 
integral equals the sum in (|E6[) for p con — > 0. The sum is 
also easily evaluated in the limit p con — > 1. In that case 
it evaluates to 



n , „ 

5>2* = ^n(n + l). 



(E9) 



Using that Li2(x) tends to — 7r 2 /6 — ln 2 (— x)/2 for large 
negative values of x [2l| , the limit of the integral 1|E7|) is 



In 2 . 2 In 27. 

n(n H . 

2 K \vl2 ' 



(E10) 



Hence for the limit of the integral to equal that of the 
sum, we require 7 = 1/V2- Inserting 7 in 1|E8|) and 
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taking the exponential, our best estimate for the product 
occurring in (|E5|) is 

n-»" 1 ( r; )■■ — ■ < En » 

j=l \ /'con; 

Using this together with £/L = 2™ we obtain (|36j) . It 
can be verified numerically that IjElip is in fact a very 



good approximation in our range of interest. Fig. [T2] 
shows a plot of the relative error as a function of n and 
Peon- For n < 45 the relative error never exceeds 3% for 
any value of 

Peon • 

For non-counting detectors, a similar derivation can be 
carried out with the recursion r\ n modified slightly since 
events where two photons reach the same detector are 
now accepted as successful connections. 
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